Laser-assisted nuclear photoeffect reexamined 
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Abstract 

The S-matrix element and the cross section of the laser-assisted nuclear photoeffect are recalculated in a 
gauge invariant manner taking into account the effect of the Coulomb field of the remainder nucleus. The 
y-photon energy dependence of the laser free cross section obtained in the plane wave and long-wavelength 
Coulomb- Volkov approximations are compared. Numerically the laser-assisted partial cross sections with 
laser photon energy 2 keV and some different polarization states of y-photon of energy 3 MeV are investi- 
gated. 
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I. INTRODUCTION 



The problem of requirement of gauge invariance of perturbation calculus in matter-field inter- 
actions in atomic physics and quantum optics was a central problem in the 60's and 70's of the 
last century, and it was finally satisfactorily clarified in the late 80's yj], M\- At this time the laser 
assisted x-ray photoeffect was widely investigated and the gauge invariant calculation of the cross 
section was also made yfl. In a recent paper |4|] the effect of intense coherent electromagnetic 
field on the nuclear photoeffect was discussed. The laser-assisted nuclear photoeffect is a process, 
which is similar to the laser-assisted x-ray photo effect (x-ray absorption). Both processes are 
bound-free transitions of charged particles (protons and electrons, respectively) that are assisted 
by an intense electromagnetic field. In both cases the initial state is strongly bound so that its 
change due to the intense laser field may be neglected. Since the correspondence between the two 
processes is straightforward, one can use the results of Q3J1 in calculating the gauge invariant cross 
section of laser assisted nuclear photoeffect making the necessary substitutions and modifications 
in the formulae. 

Therefore we re-discuss the problem of the laser-assisted nuclear photoeffect satisfying the 
requirement of gauge invariance and taking into account the effect of the Coulomb field of the 
remaining nucleus. The recoil of the remainder nucleus and the initial momentum of the y 
particle are neglected. The calculation is made in radiation (pA) gauge in long wavelength ap- 
proximation (LWA) of the electromagnetic fields. In LWA the laser field has a vector potential 
A L (t) = A [cos(o»oO e i - sm(u) t) e 2 ], that corresponds to a circularly polarized monochromatic 
wave.For the unit vectors ~e \ and ~e 2 the ~e \ ■ ~~e 2 = holds. The corresponding amplitude of 
the electric field of the laser is F = co A /c. The unit vectors of the frame of reference used are 
~e x = e\,~e y = ~e 2 and~e z = ~e 1 x~~e 2, with~e z in the direction of propagation of the intense field. 
The vector potential describing the gamma radiation is A y = Jlnhf {Vcj^Js exp [-iu> 7 ij in the 
LWA. Here tico y is the energy and s is the unit vector of state of polarization of the y photon, and 
V is the volume of normalization. The polar angles of ~~s are and <J>. The space dependent part 
of the initial nuclear (protonic) state foCr) = (2n)~ l /S 3/2 e~^ r / {fir), with fi = vy/2mE b /fi, where 
m is the rest mass of the proton and -Et, is its initial energy. The actual values of Ej, and v are 
E h = 0.137 MeV and v = 1.84 Jj. 
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II. GAUGE INVARIANT S-MATRIX ELEMENT IN THE LWA COULOMB-VOLKOV MODEL 



The wave function of a free proton in a repulsive Coulomb field of charge number Z has the 
form ip(r) = e'^xiQ^)/ W DJ. Here Q is the wave number vector of the proton. Applying 



the LWA in x(Q,r), i.e. taking X (Q, 0) 



Xc(Q), 



with 



xdQ) = 



-1/2 



(1) 



(2) 



Here a/ is the fine structure constant and A p = ^ is the reduced Compton wavelength of the 
proton. The solution (OQ) is a LWA of the Coulomb-solution. This approximation leads to the 
Fermi correction factor of the decay rate of the fi decay flg]. 

For the time dependent state of the proton in the intense field an approximate non-relativistic 
solution ifj of the time dependent Schrodinger equation of a particle in the laser plus Coulomb 
field is used, which is called Coulomb- Volkov solution D7D, [8|]. For ^ (~r , ij we use the LWA of 
nonrelativistic Coulomb-Volkov solution 



^tj,t) = V- l/2 e^xc(Q)oxp(-iEt/fl)f(t) 



(3) 



with E = fi 2 Q 2 /(2m) + U p , that is the energy of the outgoing proton in the intense field, where 
U p = e 1 F 2 l(2mtofy is the ponderomotive energy. The polar angles of the wave number vector Q 
of the outgoing proton are # and r] Q . The function fit) = exp[za s'm(oj t + rj )] where a = a & sin (■&) 
with 0!$ = eF QI (mwj). 

The gauge independent S-matrix element can be obtained with the aid of Eq.(27) of PD as 

xdQ) 



Sfi = 



vv 



J exp[i (E + E h ) t/n]f(t)jG [t (*)] ^, 



(4) 



where GLin = J (p (~r)e Iq r d 3 r is the Fourier transform of the initial stationary nuclear state, it 
is G ("7? ) = 2 yj2nfi {q 2 + ft 2 } 1 in our case, and ~~q (t) = ~Q - j^A with A = ~Al (0 + ~A T Using the 
d t G = (d q G) Eg (d A] q) (d,Aj) identity <9,G = {dfi) £ {q ■ ~E - £l ■ ~e) with £ = -±<9 ? A, i.e. 
~E = ~E L {t)+~E y . 

Now we deal with last factor of d t G. The Q ■ E L term can be neglected if the pure intense 

field induced proton stripping process is negligible since this term describes the process without 

— > — > — > — > 

the gamma photon. Furthermore, the ratio of the amplitudes of A y ■ E L and A L ■ E y equals 



iOLl^y ^ 1. Therefore the Q • E - • E = Q- E y - ^ A L - E y approximation is justified to use, 
where E y = i ^2nh,a) y /Vs exp {-ico y i}. The relative strength of the Q • E y and j- A L ■ E y terms 
is characterized by the parameter 8 = eA Q / (hcQ). In the laser free case Q = ^jlm \^ioj y - E^/fi. 
Numerical estimation shows that 6 < 0.004 if %a> y > 3 MeV and 6 increases up to 6 - 0.05 
if fia> y = 0.189 MeV in the case of laser photon energy and intensity values discussed in [Q]. 
Therefore the Q ■ E y term is the leading one in the last factor of d t G that after the substitution of 
the concrete form of d q G results 

9. \dq )t, q W+F? % 



As to the denominator of ©, the effect of A y is negligible m~q(t) and thus ~~q (t) = Q- f^A L . It 
was shown above that the amplitude of oscillation of ~q (t) due to the intense field can be neglected. 
Moreover the amplitude of oscillation of q (t) compared to fi is less than 1.4%. Therefore q = Q 
can be used in the denominator of ©. 

Using the Jacobi- Anger formula in the Fourier series expansion of f*(t) [[9Q the S -matrix ele- 
ment can be written as 

A 2ni (dG\ e JTjAoL 

Sfi = 2_j ~Y 6 MQ)]Xc(Q) h- ~ n ~ M„ (i, a) (6) 

n=no ' ' 1=Q 

with M„ (£, a) = i;J n (a)e~" m , where £ = Q ~s /Q, J n (a) is a Bessel function of the first kind, and 

JiQ 2 U p + E b 
= — + Z w r _ na) o- (V) 



III. GAUGE INVARIANT CROSS SECTION OF LASER-ASSISTED NUCLEAR PHOTOEF- 
FECT 



The cross section has the form 



where the partial cross section 



= JVn, (8) 



n=no 



Vn = 0- n{ ){Q n )\l^ n \ 2 (9) 



with 



CTnoiQn) =XciQn)a ^" 



2nA p 



(10) 

q=Q„ 



crnoiQn) = l6a f ^ c {Q n ) 7^^ . (11) 



In our case 

Here <2„ = | ^2m [ft (o» r + noj^j — U p — £/,], k y = o) 7 /c. The cases < and n > correspond to 
laser photon emission and absorption, respectively. n (< 0) is the smallest possible value of n, it 
just fulfills the h (co y + nco<^ - U p - E b > condition (\n \ ^o> 7 - t/j, - E fc J /ha>o). 



fj z n (a n )sm&d&drjo, (12) 

o Jo 

where ar„ = ar„# sin# with a n & = eF Q„f (mafy. 

The partial cross section is proportional to the Coulomb factor Xc(Q n ). The Coulomb factor 
and therefore also the partial cross section rapidly decrease with the decrease of Q n , i.e. near the 
threshold (fiajy — » Eb), and it is the fact in the laser free case too. Therefore near the threshold the 
cross section is unobservable. 

If the direction of ~s in the gamma flux is random then |/j„| 2 must be averaged as {\pi n \ 2 ) = 
J|/i„| 2 sin0J0JO/ {An). Applying the spherical harmonics addition theorem, the orthonormal 
properties and the sum rule of the spherical harmonics [10\ one can obtain i^ 2 ^ = 1/3 and 

J -/*(c^sin#)sin#^77 . (13) 

Than the average of the partial cross section cr n reads (cr„) = cr n0 (Q n ) (\fi n \ 2 y 

It can be seen from (fT3l) that the averaged differential partial cross section can be written as 



d(cr„) 1 2 . 

— = cr„ -/,,(«„# sin #), (14) 
ai2 e 3 

where cIQq = sin &d&dr] . 

In the n ^ channels the proton emission vanishes in the direction parallel {& = or & = n) 
with the laser beam since lim^o J n (x) = 0. In the n = case in this direction = <x 00 | since 
lim^o-A) (x) = 1. 

In calculating (cr n ) one can use the 7_„(x) = J„(-x) = (-1)" J n (x) relations to change the C d& 
to 2 f* /2 d&, and one can apply the 

•vr/2 

/ 2 [ar M # sin(#)] sin(#) J?? = 

ar„# 



r 

Jo 



7 2 [a m? sin {■&)] sin (#) d# = V J 2n+2k+l {2a n& ) (15) 



k=0 



[11] and the 



2 2_ J J 2n+ zk+i(2a!n»)= I J ln {x')dx' (16) 



0] 



identities. Using the x' = 2a n $x and dx' = 2a n& dx change of variable 

An f 1 

(o-„) = — cr n0 J 2n (2a n& x)dx. (17) 



In the weak field limit (a n ® — > 0) the J Q J 2n (2a n #x) dx — > J Q 7 2; , (0) = because 
lim.^o hn W = in the case of n ± and the J Q J 2n (2a n $x) dx — > J () 7 (0) = 1 because 
lim A ^ -A) (•*) = 1 in the n = case. Therefore the averaged total cross section (cr) — » ycr o. On 
the other hand (cr) —> <t t , that is the total cross section of the laser free case and with random 



polarization in the y-flux therefore cr o = ^cr r . 

If we are far from the threshold, i.e. if fioj y » Eb, then the Q„ = Q approximation holds for 
the intensity and photon energy parameter pairs (7 = 10 18 W/cm 2 with fiojQ = 2 eV; I = 6.25 x 10 21 
W/cm 2 with Hcoq = 100 and 200 eV; I = 4.0 X 10 21 , 1. x 10 23 and 2.5 x 10 24 W/cm 2 with hco = 2 
keV) discussed in |4|] and one can use the o~„o — Cqq — ~^cr j substitution in all the above formulae 
of the partial cross section. Using this approximation in (fl4l) and applying 7 2 (x) = 1, since 

n is a negative integer of large magnitude, ^ - -^cr T and (cr) ^ cr T at the intensities and laser 
photon energies discussed. 

Now it is supposed that the y-photon is polarized. If the polarization vector ~e = ~e 2 , 

the polarization vector ~s lies in the plane of polarization of the circularly polarized laser beam 

— * — > — > 
(case pol, 1) then £ = Q ■ e 2 /Q = sin # sin rj . If the polarization vector e of the y-photon is 

~~e = ~e i x~e 2 (case pol, 2) then £ = cos ft. Using (PT2l) and carrying out the integration over rj Q 
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1/4 It* i = * sin 0) sin J 0rf0, (18) 



\Mn\l i,2 = I I 7 2 (a # sin#)cos 2 #sin#d#d?7o, (19) 
Jo Jo 



and <x„ has the form 

3 

0- n ,pol,j = -^0- T \Mn\ 2 pol,j ■ ( 2 °) 

The total cross sections cr po i j = °~n,poi,j can be obtained applying £^ 7 2 (x) - 1 again resulting 
o~ P oi,\ = CTpoia = °~t valid in the cases of the laser intensities and photon energies discussed. 



IV. NUMERICAL RESULTS 



First the E p = fra) y -E b dependence of the laser free averaged total cross section is investigated. 
E p is the kinetic energy of the outgoing proton, fioj y is the energy of the y photon and E b = 0.137 
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FIG. 1 : The full line shows the E p - haj y - Eb dependence of the laser free averaged total cross section 
given by cr T = ytroo (for cr o see (fTTT ) with Qo). E p is the kinetic energy of the outgoing proton, fia) y is the 
energy of the y photon and Ey - 0.137 MeV is the binding energy of the proton initially bound in 8 B. For 
comparison the ^r-croolXc (Qo) is also plotted as a dotted line. 



MeV is the binding energy of the proton initially bound in 8 5 [|4J], Ill2ll . The charge number of the 
final nucleus ( 7 #e) is Z = 4. The full line in Fig. 1. shows cr T = ycr 00 , which is the laser free, 
averaged cross section of our model (for <t o see (fTTT) with Q ). For comparison the y cr 00 /;^ 2 (Q ) 
is also plotted as a dotted line. That is the result of the gauge independent laser-free calculation 
in the plane wave approximation, i.e. without the Coulomb correction. On the base of Fig. 1. 
one can conlude that the Coulomb correction becomes more essential with decreasing y photon 
energy. 

Next the averaged partial cross sections (<x„) (applying the cr nQ = ^cr T substitution in (fT71) ) are 
investigated numerically with fico = 2keV, hoj y = 3 MeV and the intensities discussed in [4], i.e. 
at / = 4.0 x 10 21 , 1. x 10 23 and 2.5 x 10 24 W/cm 2 (Fig. 2). Figs. 3 and 4 show the partial cross 
sections <x„ pol j in the two cases of polarization of the y photon discussed. Fig. 3 is devoted to the 
case of"e = ~e 2 (case pol, 1) and Fig. 4 shows the partial cross sections in the case of ~e =~e \ x~e 2 
(case pol, 2). 

The laser free, averaged total cross section cr r = 1.02 mb at fioj y = 3 MeV. At this y energy 
and in all cases we obtained (a) = cr po i\ = cr po / 2 = a T = 1.02 mb contrary to cr = cr n ~ 63.4 
mb obtained by |4|]. Also at this y energy the gauge independent laser- free result in the plane wave 
approximation gives <r T /Xc (Qo) = 4.11 mb. Thus one can conclude that the 62.2 times larger 
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total cross section of [4] compared to our result may be originated partly from the gauge invariant 
calculation (a factor of 15.4) and partly from ignoring the Coulomb repulsion between the final 
particles by D4J] (a factor of 4.03). In the E p — > {fno y — » 0.137 MeV^j limit this type of difference 
is more enhanced. 

V. SUMMARY 

The problem of laser-assisted nuclear photo-effect was discussed in a gauge invariant manner 
and taking into account the effect of the Coulomb repulsion between the ejected proton and the 
remainder nucleus. The investigation is mainly concerned with those y photon energies that are 
far from the threshold. In our model the calculation of the total cross section leads to 0.0161 times 
smaller result than that of [4] at the y photon energy %cj 7 = 3 MeV and in all cases of state of 
polarization of the y photon discussed. It was found that the hindering effect of the Coulomb 
repulsion in the final state, that is manifested in the appearance of the Coulomb factor in the cross 
section is huge for small kinetic energies of the outgoing proton and it must also be taken into 
account in the MeV energy range. At %co y = 3 MeV it causes about a factor of 1 /4 decrease of 
the total cross section, which is incorporated in the factor 0.0161, compared to the result of the 
plane wave approximation. The further 0.065 times decrease of the total cross section may be 
originated from its gauge invariant calculation. The numerical investigation, similar to |4|], was 
made at fraj y = 3 MeV, and it shows, that the main effect of the presence of the laser field is that 
the total cross section is distributed between the partial cross sections of the channels absorbing or 
emitting different numbers of laser photons. 
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FIG. 2: The averaged partial cross sections (<x„) for some intensity parameters discussed in Q4[| ; the inten- 
sities are (a) / = 4.0 x 10 21 (6)1. x 10 23 and (c) 2.5 x 10 24 W/cm 2 with photon energy %ojq = 2 keV. The 
polarization vector of the y photon of energy 3 MeV is random. 
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FIG. 3: The partial cross sections <t„^ / i in the case of the state of polarization of the y photon e - ei. 
The intensities are (a) / = 4.0 x 10 21 (»1- x 10 23 and (c) 2.5 x 10 24 W/cm 2 with photon energy ha> = 2 keV. 
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FIG. 4: The partial cross sections cr n po i2 in the case of the state of polarization of the y photon e = e i x e 2- 
The intensities are (a) / - 4.0 x 10 21 O)l. x 10 23 and (c) 2.5 x 10 24 W/cm 2 with photon energy hco = 2keV. 
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